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ABSTRACT 

Recent studies of lensing clusters reveal that it might be fairly common for a galaxy cluster 
that the X-ray center has an obvious offset from its gravitational center which is measured 
by strong lensing. We argue that if these offsets exist, then X-rays and lensing are indeed 
measuring different regions of a cluster, and may thus naturally result in a discrepancy in the 
measured gravitational masses by the two different methods. Here we investigate theoretically 
the dynamical effects of such lensing-X-ray offsets, and compare with observational data. We 
find that for typical values, the offset alone can give rise to a factor of two difference between 
the lensing and X-ray determined masses for the core regions of a cluster, suggesting that 
such "offset effect" may play an important role and should not be ignored in our dynamical 
measurements of clusters. 
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1 INTRODUCTION 

Galaxy clusters, the largest gravitationally-bound structures in the 
universe, are ideal cosmological tools. Accurate measurements of 
their masses provide a crucial observational constraint on cosmo- 
logical models. Several dynamical methods have been available to 
estimate cluster masses, such as (1) optical measurements of the 
velocity dispersions of cluster galaxies, (2) measurements of the 
X-ray emitting gas, and (3) gravitational lensing. Good agreements 
between these methods have been found on scales larger than clus- 
ter cores. 

However, joint measurements of lensing and X-rays often 
identify large discrepancies in the gravitational masses within the 
central regions of clusters by the two methods, and the lensing mass 
has always been found to be 2 - 4 times higher than the X-ray 
determined mass. This is the so-called "Mass Discrepancy Prob- 
lem" (Allen 1998; Wu 2000). Many plausible explanations have 
been suggested, e.g., the triaxiality of galaxy clusters (Morandi et 
al. 2010), the oversimplification of the strong lensing model for 
the central mass distributions of clusters (Bartelmann & Steinmetz 

1996) , the inappropriate application of the hydrostatic equilibrium 
hypothesis for the central regions of clusters (Wu 1994; Wu & Fang 

1997) , or the magnetic fields in clusters (Loeb & Mao 1994). 
Recently Richard et al. (2010) present a sample of 20 strong 

lensing clusters taken from the Local Cluster Substructure Survey 
(LoCuSS), among which 18 clusters have X-ray data from Chan- 
dra observations (Sanderson et al. 2009). They show that the X- 
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ray/lensing mass discrepancy is 1.3 at 3cr significance — clusters 
with larger substructure fractions show greater mass discrepancies, 
and thus greater departures from hydrostatic equilibrium. 

On the other hand, lensing observations of the bullet cluster 
1E0657-56 (Clowe et al. 2006), combined with earlier X-ray mea- 
surements (Markevitch et al. 2006), clearly indicate that the gravi- 
tational center of the cluster has an obvious offset from its baryonic 
center. Furthermore, recent studies (Shan et al. 2010) of lensing 
galaxy clusters reveal that offset between the lensing center and X- 
ray center appears to be quite common, especially for umelaxed 
clusters. Among the recent sample of 38 clusters of Shan et al. 
(2010), 45% have been found to have offsets greater than 10", and 
5 clusters even have offsets greater than 40" . Motivated by such ob- 
servations, we propose to investigate galaxy cluster models where 
the center of the dark matter (DM) halo does not coincide with the 
center of the X-ray gas (See Figure 1). 

If the X-ray center of a cluster has an offset from its lens- 
ing (gravitational) center, then the X-rays and lensing are indeed 
measuring different regions of the cluster. Given the same radius, 
the lensing is measuring the DM halo centered at the gravitational 
center (shown by the dark blue sphere in Figure 1 ), while the X- 
rays are measuring the sphere of the halo that is offset from the 
true gravitational center (shown by the red circle in Figure 1). In 
this case, there will always be a natural discrepancy between the 
lensing and X-ray measured masses — or specifically, the X-ray 
mass will always be lower than the lensing mass, just as the long- 
standing "mass discrepancy problem" has indicated. 

In this paper, we investigate the lensing-X-ray mass discrep- 
ancy caused by the offsets between DM and X-ray gas. To check 
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Figure 1. Oifset between the dark matter center and the X-ray center in a 
galaxy cluster 

our predictions, we compile a sample of 27 clusters with good lens- 
ing and X-ray measurements. We conclude that such "offset" effect 
should not be ignored in our dynamical measurements of galaxy 
clusters. A flat ACDM cosmology is assumed throughout this pa- 
per, where f2,„=0.3, Q.f^=Q.l, and Hq = 70kms"'Mpc"'. 



2 MASS DISCREPANCY AS A RESULT OF THE DARK 
MATTER-BARYON OFFSET 

We model our galaxy cluster with a fiducial model as the following: 
(1) the DM halo is modeled by the Navarro-Frenk- White (NFW) 
profile (Navarro et al. 1997) with concentration c = 4.32 and scaled 
radius r, = 516 kpc, (2) the gas distribution is modeled by a yS 
model with /? = 0.65, the cluster core radius = 200 kpc, and 
the gas fraction f„^^ = 12%, (3) the mass density of the BCG is 
described by a Singular Isothermal Sphere (SIS) with a velocity 
dispersion of 300km/s. 

The projected mass within a sphere of radius is 

J nil! pR^ 
[5;NFw(/^') + 2=as(^) 
Jo 

+UcGiR')\R' dR' de, 

where R' = Vrf2 +R^ + 2dR cos 6 is the 2-D radius from the halo 
center, R is the 2-D radius from the X-ray gas center, d is the 2-D 
offset between the halo center and X-ray center, and Enfw, Sgas, 
and Ebcg are the projected mass densities of the DM halo, the 
gas and the BCG, respectively. For a given radius R^, the gravi- 
tational mass measured by lensing mie„s can be given by m(/? j, 0) 
(as shown by the dark blue sphere in Figure 1), while the projected 
mass measured by X-rays njxray is described by m(/?j-, d) (the mass 
within the red circle in Figure I). We now calculate the mass ratio 
m(R^,d)/m{R^,0), or equivalently, raie„s/'Wxmy 

Figure 2 shows the mass ratio as a function of the 2-D offset d, 
for a typical rich cluster. The solid curves are the mass ratio with the 
fiducial model, the dashed and dotted curves are the mass ratio with 
the NFW concentration c = 4.04 and 5.13 (top left), the cluster core 
radius = 150 kpc and 400 kpc (top right), the p index p = 0.6 and 
0.9 (bottom left), the gas fraction /gas = 0.1 and 0.2, respectively. 
For these cases, the three curves from top to bottom are for the three 
measuring radii R^ = 50 kpc, 100 kpc, 200 kpc, respectively. From 
Figure 2 we have the following conclusions: 
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Figure 2. Ratio of projected gravitational masses, as a function of the 2-D 
offset d and the measuring radius R^-. The solid curves are the mass ratio 
for the fiducial model with c = 4.32, r, = 516kpc, /3 = 0.65, = 200kpc, 
and /gas = 0.12. The dashed and dotted curves are for the NFW concentra- 
tion c = 4.04 and 5.13 (top left), the cluster core radius = 150 kpc and 
400 kpc (top right), the /3 index /3 = 0.6 and 0.9 (bottom left), the gas frac- 
tion /gas = 0. 1 and 0.2, respectively. The three dotted (dashed, solid as well) 
curves from top to bottom in one panel con'espond to Ry = 50, 100, 200 kpc, 
respectively. 



(1) The lensing measured mass miens is always higher than the 
X-ray measured mass /Wxray For typical values of offset d = 100 kpc 
and Ry = 100 kpc, raie„s/mxray ~ 2, comparable to the ratio found in 
early studies (Allen 1998; Wu 2000; Richard et al. 2010). 

(2) The "Offset Effect" we are reporting here should contribute 
significantly to the long-standing "Mass Discrepancy Problem". 

(3) The ratio of m]^„Jm„^y increases with offset d. 

(4) miens //nxiay depends very strongly on . Here R^ acts like 
the arc radius rare in strong lensing, i.e., we only measure the en- 
closed mass within a small region of < R^ . When R^ is very small, 
the ofi'set effect is most prominent and gives large raiens/^Jxiay In- 
creasing Ry will reduce miens /"Jxiay When is very large (com- 
pared with d), the offset effect will be "smeared out", and the ;nie„s- 
"ixray discrepancy introduced by the offset will vanish. 

(5) The mass ratio is very sensitive to the NFW concentration, 
and it increases dramatically with c. 

(6) The mass ratio increases with the core radius, and de- 
creases with y0 index and gas fraction. However, the mass ratio is 
not very sensitive to the gas model. 



3 COMPARISON WITH OBSERVATIONAL DATA 

To compare with our theoretical predictions, we compile a sample 
of 27 clusters with 48 arc-like images, which have both strong lens- 
ing and X-ray measurements. The clusters and their lensing and X- 
ray data are listed in Table 1. For the 22 arcs that have no redshift 
information, we estimate their lensing masses mie„s by assuming 
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the mean redshifts of {zd) = 0.8 and 2.0, respectively. The X-ray 
data are taken from Tucker et al. (1998), Wu (2000), Bonamente 
et al. (2006), and references therein. The offsets between lensing 
and X-ray centers are taken from Shan et al. (2010). The clusters in 
our table are classified as relaxed (with cooling flow) and unrelaxed 
(which are dynamically unmature), from their X-ray morphologies. 
The definition has been used in the literature by Allen (1998), Wu 
(2000), Baldi et al. (2007), and Dunn & Fabian (2008). 

Mass from strong lensing. Assuming a spherical matter distri- 
bution, one can calculate the gravitational mass of a galaxy cluster 
projected within a radius of Tjuc on the cluster plane as 



m|cns(< ''arc) = J^r^^^^n 



(1) 



-crit - 4;^-oj[)j- is the critical surface mass density, £);, D, 



where Ecri, = 

and D/j are the angular diameter distances to the cluster, to the 
background galaxy, and from the cluster to the galaxy, respectively. 
The above equation is actually the lensing equation for a cluster 
lens of spherical mass distribution with a negligible small align- 
ment parameter for the distant galaxy within r.^^. The values of 
'"lens within the arc radius r^y^ are listed in Table 1 . 

Allen (1998) pointed out that the use of more realistic, ellip- 
tical mass models can reduce the masses within the arc radii by 
up to 40%, though a value of 20% is more typical. However, such 
corrections are still not very significant compared with the large 
discrepancies between the lensing and X-ray determined masses. 
We will discuss it in more detail in the next section. 

Mass from X-rays. Assuming that the intra-cluster gas is 
isothermal and in hydrostatic equilibrium, the cluster mass m(r) 
enclosed within a radius r can be easily calculated from 



Gm(r) kT dln«gas('-) 



fim,, 



dr 



(2) 



where T is the gas temperature, «gas the gas number density, the 
proton mass, and yu = 0.585 the mean molecular weight. Here we 
assume that the gas follows the conventional /? model, i.e., ngas('') = 
"Eas(0)(l -I- r-/r^) T . In order to compare the mass measured by X- 
rays with the lensing result, we need to convert this m(r) (i.e., 3-D) 
into the projected mass m„^y (see e.g. Wu 1994): 



where 



1.13 X lO'W 



f(y) = 



R\l r, M kT 



'-J\0.1Mpc/\lkeV 



(3) 



(4) 



2(1 +y2)l/2' 

the mass ratio "iiens/'^xray are listed in Table 1. 

Figure 3 shows the relation between the mass ratios miens /"^xray 
and the (scaled) offsets for our sample of 27 clusters (48 arc im- 
ages). It should be pointed out that the 27 clusters in our sample 
have quite different sizes and masses. This can be seen from the 
wide range of the cluster temperatures — from 4 keV to 14 keV. 
Therefore, it is useful to compare the offsets of the clusters on 
the same scale. We realize that the M-T relation of clusters scales 
asM ~ T^/^ (e.g., Nevalainen et al. 2000; Xu et al. 2001), and that 
M ~ R^, where R is the size of the cluster. Therefore, in Figure 3, 
instead of using the physical offset d, we use a scaled offset which 
is characterized by \^cl'^ki:\^^^ ■ 

From Figure 3, the mass ratios miens/"^xray exhibits large dis- 
persions — roughly ranging from 2 to 4. Many clusters have large 
error bars. It appears that relax clusters (marked by crosses) have 
smaller 'Miens/jJixray ratios. The fact that mjens > "'xiay is consistent 
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Figure 3. The ratio of lensing and X-ray determined masses for our sample 
of 27 clusters (48 arc images). The x-label show the scaled offset between 
DM and baryons. The squares denote unrelaxed clusters, and crosses the 
relaxed clusters. The dashed line shows a fit satisfying 'niens/"ixray ~ 
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0.75 for the clusters with offset larger than 



is also roughly consistent with our predictions as plotted in Fig- 
ure 2. 

However, no strong correlation has been found between the 
offset and mass discrepancies. We notice that many clusters in the 
sample have very small offset values — smaller than the errors in 
lensing and X-ray measurements which are typically a few arcsec- 
onds. So these offset values are not robustly measured themselves, 
and we thus remove these data points and only focus on clusters 
with large offsets of d > 10", as has been suggested in Shan et 
al. (2010). This leaves a sub-sample of only 24 arc images. The 
dashed line in Figure 3 shows a.x^ fit to this sub-sample, which sat- 
isfies mi^,,Jm„.^y 



3.24('-^^^f^° with a reasonable x" = 0.75. 



with our theoretical predictions, and the ratio of niiens/'^xi 
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VA-r/keV 

We can find mig„slm_„ay increasing slightly with d. 



4 DISCUSSION AND CONCLUSIONS 

As has been reported by Shan et al. (2010), it might be fairly com- 
mon in galaxy clusters that the X-ray center has an obvious offset 
from the gravitational center. We have explored the dynamical con- 
sequences of this lensing-X-ray offset and tried to attribute such an 
effect to the long-standing "Mass Discrepancy Problem" in galaxy 
clusters. Our theoretical model predicts that such an offset effect 
will always result in a larger ozics than 'Jixray, with a typical mass 
ratio JKiens/nJxray ~ 2, wMch is consistent with observations. 

To test our model, we have compiled a sample of 27 clus- 
ters, and studied in detail their lensing and X-ray properties and 
obtained their lensing and X-ray masses, mie„s and oj^ay The lack 
of strong correlation between m]^„Jm^y.^y and the offset d suggests 
that the problem is more complicated. As we have found in Sec- 
tion 2, "ijens/'Tixray is not Only a function of d, but also depends very 
strongly on Ry (or the arc radius r.^y^). Apparently, each cluster in 
our sample has quite different r.^y^. 
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Probably, other mechanisms than the offset effect should play 
important roles, and the lensing-X-ray mass discrepancy may not 
be just from one mechanism, but a combination of many effects: 

(1) The central regions of clusters may be still undergoing dy- 
namical relaxation, and the X-ray gas may not be in good hydro- 
static equilibrium. Therefore, large errors could be induced in the 
X-ray measurement of cluster cores, especially for unrelaxed clus- 
ters. 

(2) The spherical models are too simple to reflect the real mass 
distribution of clusters. The use of more realistic mass model could 
reduce the lens mass within the arc radius by up to 40%, though 
values of ~ 20% are more typical (Bartelmann 1995; Allen 1998). 

(3) The presence of substructures may complicate our simple 
spherical lens model, and hence could be a main source of uncer- 
tainties in ;niens. The absence of the secondary arc-like images in 
most arc-cluster systems may indicate the limitations of the spher- 
ical mass distribution in the central regions of clusters. 

It should be noted that the mass ratios we obtained here are 
slightly higher than Allen (1998) and Wu (2000) because they un- 
fortunately used a Hubble constant of Ho = 50kms 'Mpc The 
use of Ho = 70kms"'Mpc"' here will of course make the mass 
discrepancy problem more pronounced. 

It should be noted that the gas represents only a 10% pertur- 
bation due to the small ratio of gas-to-DM in the central region, 
likewise the offset of the gas is only a small perturbation (less than 
10%) to the otherwise concentric matter density or potential. It is 
unlikely to create a factor of two difference in the lensing-derived 
enclosed masses within an arc. 

To illustrate the lensing effect of the offset perturbation and 
triaxiality, we show the critical curves in Figure 4. The solid curves 
indicate the critical curve of circular NFW plus f} model with- 
out offset, the dotted curves indicate the critical curve of ellip- 
tical NFW plus f} model with offset d = 10". The square and 
cross denote the center of dark matter and the hot gas, respec- 
tively. For the NFW profile, c = 4.3, r, = 516kpc; for the fi model, 
P = 0.65, rc = 150kpc. We also introduce the triaxiality with the 
ellipticity e = 0.15 and position angle 6 = 30°. We also assume 
the lens and source redshifts zi = 0.3, z, = 1. We can see that the 
predicted critical curves (dotted lines) have very similar sizes as the 
predicted critical curves for a benchmark model (solid lines) with 
the same mass DM and gas mass but in concentric spheres. 

Early studies have suggested that statistically unrelaxed clus- 
ters have larger mass discrepancies than relaxed clusters (Allen 
1998; Wu 2000; Richard et al. 2010). As Shan et al. (2010) have 
reported, the clusters with large offset of d > 10" are all unrelaxed 
clusters. If such offsets exist and are big, then they must come into 
play in our dynamical studies of galaxy cluster, and should not be 
ignored, especially for unrelaxed clusters. 
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Table 1. The X-ray and lensing mass discrepancies of 27 clusters. For the 22 arcs that have no redshift information, we estimate the mean redshifts of {zd) = 0.8 
and 2.0, respectively. Refs A and B give the references of the lensing and X-ray data, respectively. The last column shows the classification of the clusters: 
"R/U" means relaxed/um'elaxed. 
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7 on+0.68 
^■°-^-0.61 


R 


13 


A 1689 
A1835 


0.183 
0.252 


0.60 
1.61 


1.85 
6.33 




0.20 
0.17 


4.5(3.8)'' 
2.82(2.23)'' 


8 
11 


n (11 +0.40 
-0.30 

9 8+' '' 
^■"-1.4 


0.65 +".''4 
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"■"■'-0.04 


14+°.°2 

o.o8:J||{ 


i.68;»:« 

1 79-1-0.36 
-0.36 


2 58+038,2 29+°32) 

1 70+°.'*(l 35+°.2*) 
-0.36' -0.28' 


R 
R 


13 
13 


A1914 


0.171 


11.3 


32.9 




0.10 


1.16(1.01)'' 


10 


a g+0.3 
-Q.3 


90+°'" 
-Q.OL 


0.200!»;°| 
"■^'^-0.0003 


Q 70+0.036 
^'■'^-0.036 


157+0.09,1 44+0.07) 

-P.99 -P.97 
^^l+0.O3.^^;^+0.O2^ 

"■^1-0.03 ("■''3-0.02' 
O49O.02 ,0 43+0.02) 
-0.02 ' -0.02' 


U 


2 


A2204" 


0.151 


1.20 


3.15 




0.025 
0.01 


0.08(0.07)'' 
0.013(0.012)'' 


10 




48+™- 
''■*°-0.002 


p. 1 1 -1-0.0040 
^■^ -0.0040 

n n^+00009 

-p.0QQ9 


R 


2 


A2163 


0.203 


44.0 


146.9 


0.73 


0.07 


0.58 


1 


14 6*°-^^ 

7 1+°-- 
'■'-0.2 


62+°°- 
"■''^-0.02 


33+0-'*2 
"■^^-0.02 


Q 23+0.033 
-0.033 


2 30+0.3; 


U 


13 


A2218" 


0.176 


19.1 


56.9 


1.03 
0.70 


0.28 
0.09 


8.60 
0.89 


11 


n ^r+0.08 
0.65 


25"^°-^^ 
^■^■^-0.05 


1.67_„3, 
0.25+°." 


5 07+'^''' 
J."'_0.93 
3 94+1.73 
^■'+-1.04 


U 
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2.52 


0.09 


0.82 










n T^+0 11 

"■25-0.065 

0.38+°3| 


1 -)n+l 40 
^■2"-0.S5 

2 19+1.17,1 78+0.95) 
-' -l.n''^'°-Q.95' 
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0.228 


11.3 


41.2 
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0.79(0.64)'' 


11 
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40+°°' 
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0.164 
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0.04 


0.13 


10 


s ^,+0.3 


58+°.°- 
"■=°-0.02 


"■'4tSl 
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U 


2 
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1.31 
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0.12 
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10 
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0.228 


6.00 


21.9 


0.91 


0.20 


3.8 


10 


11 1 + '° 
"■'-1.0 


59+°°2 
"■^^-0.02 
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^■^^-0.32 


U 


13 


CL0024 


0.395 


13.2 


70.4 


1.68 


0.26 


4.7 


6 
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48+°°" 
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U 


13 
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1.50 
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1 
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13 
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2.79 
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10 
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■ -0.02 
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R 


13 


MS 1455 
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2.77 


11.1 
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1.22(0.96)'' 
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^■^^-0.28 
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07+°°' 


5 7+0.077 
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U 


13 


MS2053 
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3.15 
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10 
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"■' -0.01 
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2 47+0.54 


U 


2 
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5.70 


26.1 
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0.99(0.72)'' 
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"■"■'-O.OI 
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-0.48''^''^'-a35' 


R 


13 
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0.42 


10 
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06+°°' 
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3 07+0.54 
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R 


13 
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0.451 
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0.28 
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1 
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R 


13 



" Multiple-arc system. 
Arc-like image is assumed at Zs = 0.8 (zs = 2). 
R: Relaxed, U: Unrelaxed. 

References: (I) Allen 1998; (2j Bonamente et al. 2006; (3) Bradac et al. 2006; (4) Clowe et al. 2006; (5) Gioia et al. 1998; (6) Jee et al. 2007; (7) Kneib et al. 1993; (8) Limousin et al. 2007; (9) Newbury & 
Fahlman 1999; (10) Sand et al. 2005; (1 1) Smith et al. 2005; ( 12) Tucker et al. 1998; (13) Wu 2000 



